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Abstract 

Background: The accumulation of deleterious mutations of a population directly contributes to the fate as to 
how long the population would exist, a process often described as Muller's ratchet with the absorbing phenomenon. 
The key to understand this absorbing phenomenon is to characterize the decaying time of the fittest class of the 
population. Adaptive landscape introduced by Wright, a re-emerging powerful concept in systems biology, is used 
as a tool to describe biological processes. To our knowledge, the dynamical behaviors for Muller's ratchet over 
the full parameter regimes are not studied from the point of the adaptive landscape. And the characterization of 
the absorbing phenomenon is not yet quantitatively obtained without extraneous assumptions as well. 

Methods: We describe how Muller's ratchet can be mapped to the classical Wright-Fisher process in both 
discrete and continuous manners. Furthermore, we construct the adaptive landscape for the system analytically 
from the general diffusion equation. The constructed adaptive landscape is independent of the existence and 
normalization of the stationary distribution. We derive the formula of the single click time in finite and infinite 
potential barrier for all parameters regimes by mean first passage time. 

Results: We describe the dynamical behavior of the population exposed to Muller's ratchet in all parameters 
regimes by adaptive landscape. The adaptive landscape has rich structures such as finite and infinite potential, real 
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and imaginary fixed points. We give the formula about the single click time with finite and infinite potential. And 
we find the single click time increases with selection rates and population size increasing, decreases with mutation 
rates increasing. These results provide a new understanding of infinite potential. We analytically demonstrate the 
adaptive and unadaptive states for the whole parameters regimes. Interesting issues about the parameters regions 
with the imaginary fixed points is demonstrated. Most importantly, we find that the absorbing phenomenon is 
characterized by the adaptive landscape and the single click time without any extraneous assumptions. These 
results suggest a graphical and quantitative framework to study the absorbing phenomenon. 



Background 

Muller's ratchet proposed in 1964 is that the genome of an asexual population accumulates deleterious 
mutations in an irreversible manner. It is a mechanism that has been suggested as an explanation for the 
evolution of sex [l] . For asexually reproducing population, without recombination, chromosomes are directly 
passed down to offsprings. As a consequence, the deleterious mutations accumulate so that the fittest class 
loses. For sexually reproducing population, because of the existence of recombination between parental 
genomes, a parent carrying high mutational loads can have offspring with fewer deleterious mutations. 
The high cost of sexual reproduction is thus offset by the benefits of inhibiting the ratchet pi. Muller's 
ratchet has received growing attention recently. Most studies of Muller's ratchet are related to two issues. 
One is that without recombination, the genetic uniformity of the offspring leads to much lower genetic 
diversity, which is likely to make it more difficult to adapt [3]. So its adaptiveness arouses concern. The 
other is that population lacking genetic repair should decay with time, due to successive loss of the fittest 
individuals [4] [5]- So the fixation probability arouses concern. In addition, Muller's ratchet is relevant 
to some replicators [f| [7j, endosymbionts |8j, and mitochondria [9j. In order to assess the relevance of 
Muller's ratchet, it is necessary to determine the rate (or the time) for the accumulation of deleterious 
mutations [10] . It is widely recognized that the rate of deleterious mutations being much higher than that 
of either reverse or beneficial mutations results in a serious threat to the survival of populations at the 
molecular level [21. Because models proposed must rest on the biological reality, which must be analyzed 



on their own without any injection of extraneous assumptions during the analysis 11 . Overall, it has been 
a long interest to develop a suitable and quantitative theory for the ratchet mechanism and the incidental 
absorbing phenomenon. 



Biologists have suggested 12 10 that a quantitative framework is needed. The potential evolutionary 
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importance of Muller's ratchet makes it desirable to carry out careful quantitative studies 12 . And the 
incidental absorbing phenomenon is investigated quantitatively in broad literature. The simplest and earliest 
mathematical model is the pioneering work in Ref. [13] . It described the same evolutionary process on 
the condition of deterministic mutation-selection balance according to the Wright-Fisher dynamics. And 
it indicated numerical evidence of relation between the total number of individuals and the average time 
between clicks of the ratchet, but it did not focus on the absorbing phenomenon. It treated the pioneering 
model as a diffusion approximation [14] , and produced more accurate predictions over the relatively slow 
regime. It noted that the increasing importance of selection coefficients for the rate of the ratchet for 
increasing values of the total number of individuals. But it is represented as stochastic differential equations 
and did not get the predictions over all parameters regions. It employed simulation approaches to Muller's 



ratchet 15 and estimated how different between the distribution of mutations within a population and a 
Poisson distribution. But it did not emphasize the absorbing phenomenon. In Ref. [2] it obtained diffusion 
approximations for three different parameter regimes, depending on the speed of the ratchet. The model 
shed new light on [14] . But it mainly focused on the property of the solution for these stochastic differential 



equations. In Ref. 10 it mapped Muller's ratchet to Wright-Fisher process, and got the prediction of the rate 
of accumulation of deleterious mutations when parameters lie in the fast and slow regimes of the operation 
of the ratchet. But it put the constraints of Dirac function on the boundary. Previous works mainly focused 
on the parameter regimes with lower or higher mutation rates. And models are represented as stochastic 



differential equation. In Ref. 16 authors imagined the population evolved on an adaptive landscape, but 
they could not analytically construct it. It described discrete birth-death model and its corresponding 
diffusion manner by the adaptive landscape [17] . But it did not discuss the absorbing phenomenon. The 
concept of adaptive landscape is proposed by Sewall Wright to build intuition for the complex biological 



phenomena [5]. In the present article, inspired by 16 10 , we model Muller's ratchet as a Wright-Fisher 
process, analytically construct the adaptive landscape, where the non-normalizable stationary distribution 
occurs. Here the adaptive landscape is analytically quantified as a potential function from the physical point 



of view 18 . We give the adaptive and unadaptive states for the whole parameters region by the adaptive 
landscape. We give the formula for the single click time of Muller's ratchet in the face of infinite and finite 
potential. In addition, we can handle the absorbing phenomenon without extraneous assumptions. 

The key concept in constructing the adaptive landscape is of potential function as a scalar function. 
There is a long history of definition, interpretation, and generalization of the potential. Such potential has 
also been applied to biological systems in various ways. The usefulness of a potential reemerges in the current 
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study of dynamics of gene regulatory networks 19 , such as its application in genetic switch 20 21 22 23 . 
The role of potential is the same as that of adaptive landscape. In this article, we do not distinct them. 

We now make the obvious advance to Muller's ratchet. We analytically construct the adaptive landscape. 
We demonstrate the position and adaptiveness of fixed points. This makes the dynamical behaviors of the 
population to be investigated. In addition, we give the area with imaginary fixed points. This makes the 
explaining for the imaginary fixed points biologically possible. Infinite potential barriers can be crossed over 
under some cases. We handle the absorbing phenomenon without any extraneous assumptions under the 
condition of diffusion approximation. Inversely, we demonstrate the power of the adaptive landscape. 



Methods 

Discrete model and absorbing boundary 

We consider here in population genetics an important and widely applied mechanism- Muller's ratchet. 
It is the process by which the genome of an asexual population accumulates deleterious mutations in an 



irreversible manner 24 25 . It corresponds to the repeated irreversible loss of the fittest class of individuals 



because of the accumulation of the deleterious mutations, the effective absence of beneficial mutations, 



without any recombination, but with the random drift 25 12 . Consider a population of haploid asexual 
individuals with discrete generations t — 0, 1, 2, . . .. The common point in a generation is regarded as an 
adult stage, after all selection has occurred and immediately prior to reproduction. New mutations occur at 
reproduction and all mutations are assumed to deleteriously affect viability but have no effect on fertility. 
Supposed population size is fixed for each generation. The viability of a newly born individual is taken to 
be determined solely by the alleles they carry. This allows us to divide the population into different classes 
with different genotypes. 

Here in one dimensional case, we consider one locus with two alleles (for example, A and a), that is, 
there are two classes in the haploid asexual population, one class with allele A while the other with allele 
o, supposed mutation from allele A to a is deleterious. We assume fixed population size of N , which means 
we have N alleles in all. We also assume that N > 1. Generations are non-overlapping. The lifecycle of 
the individuals in the population is from adults to juveniles, during which we consider irreversible mutation, 
selection, and random drift. The frequency of the allele A for generation t is p while that of allele a is 1—p. 
Let fi be the probability that an offspring of an adult with allele A is an individual with allele a, labeled by 
Mi o, that is M 10 = fj,. Analogously, M = 1 — fi, M 01 = 0, Mi : i = 1. The relative viability of individuals 
with allele A is vq = 1 while that of individuals with allele a is v\ = 1 — a. Where a can be treated as an 
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effective selection coefficient associated with deleterious mutations. So the values of parameters for /j, and a 
are from to 1. Then in generation t+1, when selection and deleterious mutation are active, the probability 
that the offspring of a parent with allele A is chosen to be with allele a is fip(l — er), the probability that 
the offspring with allele A is (1 — /j,)p, the probability that the offspring of a parent with allele a is still with 
allele a is (1 — <r)(l — p). So the frequency of allele A in generation t + 1 is 

P(t+I)= l-<T + <7(1-/X)p(i)- (1) 

Eq.(l) describes the deterministic process that ignores random drift. Under the mutation-selection balance, 
the fixed points are 

Pi(t = oo) = 0, (2) 
= oo) = (3) 

This means the population ultimately arrives at the state with allele frequency p l or f5 2 and no transition 
between the two states occur. It is evident that selection rates a is greater than mutation rates [i. But pop- 
ulations always evolve randomly. Since each individual is assigned a parent independently, if by generation 
t the average value of allele frequency is p„ )t , we have that the transition probability from p n , t to p n ,t+i by 
generation t + 1 is 

Wnm = m\(N- m)\ (P "' t)TO(1 ~ P^" 1 - n,m = 0,l,2,...,N. (4) 
The matrix of transition probabilities is W which is composed of elements W nm . The dynamical rule can 
be written as the matrix Eq.(5) 

P(i + l) = WP(t). (5) 

Where P(i) represents the probability distribution of allele A. It is composed of N + 1 elements p^(t), 
n = 0,1, ... ,N, where p^(t) denotes the probability that allele A has the frequency n/N in generation t in 
the presence of deleterious mutation, selection and random drift. The matrix of transition probabilities is 

w=(;; T ). 

Here v T means the transpose of vector v. Then Eq.(5) can be expressed as the following 

{ p(t + l) )-\0 w ) \ p(t) j' ^ 
Where is a column vector with N vanishing components, and v is a column vector with N non-zero 
elements given by v m — W 0m , where m is from 1 to N. The column vector w is an N x N matrix with 
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elements w nm = W nm , where n and m run from 1 to N. By spectrum decomposition we can derive the 

maximum eigenvalue for w, the corresponding eigenvector is the quasi-stationary probability density of allele 
A. The following is its derivation process. 

p(t+l)=wp(t), (8) 

Po(* + l) =Po(t)+v T p(t), (9) 

Eq.(8) has the solution 

p(t) = w*p(0). (10) 

Because the leading large time behavior is determined by the eigenvalue of matrix w. The issue is transformed 
to solve the leading large eigenvalue of w and corresponding eigenvector. But 

N 

Po(t) = l-^fti(t) 

71=1 

= i-i t pW 

= l-I T w*p(0). (11) 

Where I T = (1, . . . , 1), there are N elements 1 of this vector. If we denote the leading large eigenvalue of w 
is Ai, 

1-poW = I T w*p(0) 

cx A*. (12) 



We can describe Eq.(8) as 



where 



p(t)w[l-pb(t)]q(t), (13) 



Then the probability density is derived by 



r wt P(°) 

q= hm -= — . (15) 
t^oo I T w T p(0) 

Where Ai is determined by 1— Po(t). But l — p (t+l) w (1 — v T q)[l— po (t)]- In the end we get Ai = 1 — v T q. 
Populations evolution is a natural and random process. We model the process as the discrete form. The 
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boundary of the discrete model is determined by its transition matrix. Generally, The transition probabilities 
from generation t to generation t + 1 for the allele frequency are expressed as 

N - ( (l-A0Pn,t-i \ m fi (l-M)Pn,t-i 



m\(N - to)! V 1 - ° + f(l - fi)p n , t -ij V 1 - (J + cr(l - /i)p„,t-i/ 
From the expression of transition probability, it can be seen that the transition probabilities are zero for 
any frequency p n ,t-i under the condition of parameter a = 1. It means the population stays at its initial 
state. In addition, the transition probabilities from the boundary to its next are (1,0,..., 0) T . This means 
boundary can not output any probability flow to its next, it only absorbs probability from next. We call 
absorbing phenomenon occurring at the boundary 0. 

Continuous model and adaptive landscape 

Diffusion approximation 

Here we briefly outline the diffusion approximation from the discrete to continuous models. At generation 
t the frequency of allele A is i/N, after evolutionary force, at generation t + 1 the allele frequency becomes 
j/N. Here St — 1, the probability that allele frequency becomes j/N is 

N 

pj(*+i) = E^(*)- 

i=0 

Where Wij is the transition probability. Diffusion approximation is a description of the process, valid when 
N is large, where the allele frequency n/N are replaced by real number x, < x < 1. Given that A starts 
out at gene frequency Xq. Let p(x — Sx, t) be the probability for allele frequency A after t generations. And 
p(x, t + 1) be the probability of allele A after t + 1 generations, then 

p(x, t + 6t) = ^2 W ( x ' 1 + St \ x - fa > l )Pi x - Sx > l )- ( 18 ) 

8x 

Letting p(x, t) — p(x, t)/N, among this p(x, t) is the probability density. Define M(x) as the probability that 
x increases by systematic force that include mutation and selection. And define V(x) as the probability that 
x changes because of random drift, either decreasing by amount Sx with the probability V(x)/2 or increasing 
by the amount Sx with the probability V(x)/2. 

1 f 1 

M(x,t) = lim — / SxW(x,t + St\x-Sx,t)d(Sx), (19) 
st^o St J 

i r 1 

V(x,t) = lim- / (Sx) 2 W(x,t + St\x-Sx,t)d(Sx). (20) 
(St— s-o St J 
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In any time interval St, the probability that x remains at x is 1 — M(x) — V(x). The changes in states are 
only to Sx or — 5x. So Sx is 0, positive or negative. 

p(x,t+8t) = (l-M(x)-V(x))p(x,t)+M(x~5x)p(x-Sx,t) + ^V(x+5x)p(x + Sx,t) + ^V(x-5x)p(x-Sx,t), 

(21) 

and we can treat Eq.(21) as the following 

p(x,t + St) — p{x,t) = —[M(x)p(x,t) — M{x — Sx)p(x — Sx,t)} 

+~[V(x + Sx)p(x + Sx,t) - V(x)p(x,t) 
-V(x)p(x, t) + V(x - Sx)p(x - Sx, t)], (22) 

that is 



p(x,t + St)~ p(x,t) _ 6[M(x)p(x,t)} 1 5(5[V(x)p(x,t)]) 
St ~ ~ S~x + 2~ 5{5x) ' 



(23) 



Because the function about the change of allele frequency in one generation is continuous and smooth enough, 
under the condition that |pj,t+i — Pi.t\ > 1/N is of order of magnitude small than 0(1) in one generation. 
Put it differently, the change of allele frequency is not more than 1/N, and the probability density is smooth 
enough during the time scale of one generation. So we represent the process as the following approximated 
diffusion equation 

= ~[M{x)p{x,t)] + ~~[V(x)p(x,t)}, (24) 
and according to the definition of M(x) and V(x), the explicit expressions of them are 

(1 - p)x 



M(x) 

_[_ — U ~T U [ A. — j_L J Jli 

tM rr — i A — rr ( 1 — // Wl 

(25) 



1 — a + <r(l — p)x 
x[(a — p) — cr(l — p)x] 
1 — a + er(l — p)x 



V{x) = X ^ w ^. (26) 



Among this M(x) is the symbol for the change in allele frequency 26 11 that occurs in one generation due 



to systematic force. The function V(x) is the variance in allele frequency after one generation of binomial 



sampling of N alleles 27 . 



Adaptive landscape 

Under the general diffusion approximation, frequency p n t is treated as continuous quantities x, and this 
leads to the distribution of the frequency for the allele A being the probability density. Let p(x,t) be the 
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probability density of the frequency for the allele A being x at time t. The diffusion process can be expressed 
by the following symmetric equation 

8tp(x, t) = 8 X [eD{x)d x - f(x)]p(x, t) (27) 

with 

fix) = M(x)-eD'(x), 
2eD{x) = V{x). (28) 

With a prime denoting differentiation of a function with respect to its argument such as D'(x) = d x D(x). 
Where M(x) and V(x) is from Eqs.(25) and (26) respectively. Adaptive landscape is directly given when we 
consider natural boundary as Feller's classification. It is 

^=Iw)^ (29) 

The symmetric Eq.(27) has two advantages. On the one hand, the adaptive landscape is directly read out 
when the detailed balance is satisfied. On the other hand, the constructive method is dynamical, independent 
of existence and normalization of stationary distribution. We call f(x) directional transition rate, integrating 
the effects of M{x) and the derivative of V(x). Directional transition rate can give equilibrium states when 
it is linear form. 

When the process lies at stationary state, the probability flux of the system is zero, and probability flux 
flows in x G [0, 1]. In general, the stationary distribution for the diffusion approximation satisfying natural 
boundary condition is given by 

I /$(x) 

p{X,t = OOj = 7F CX P|— — 

It has the form of Boltzmman-Gibbs distribution [28], so the scalar function $(x) naturally acquires the 
meaning of potential energy [19) . The value of Z determines the normalization of p(x,t = oo) from the 
perspective of probability, and the finite value of Z manifests the normalization of p(x,t = oo). The 
stationary distribution is not true in the face of infinite Z . It demonstrates absorbing phenomenon occurs 
at the boundary. Together with the flux at the boundary, the true stationary distribution could be got. The 
constant e holds the same position as temperature of Boltzmman-Gibbs distribution in statistical mechanics. 
But it does not hold the nature of temperature in Boltzmman-Gibbs distribution. 
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We are interested in the dynamical property of adaptive landscape, so we treat $ and $/eno difference in 
this respect, that is, for convenience we can take e = 1 of eD(x). So according to Eq.(29), we have adaptive 
landscape as the following 



From the expression of adaptive landscape $(x), we may find there are two singular points and 1 of 
adaptive landscape, characterized by infinite value, infinity means adaptiveness or unadaptiveness of the 
system. Here the adaptive landscape is composed of three terms. The first term and the third term quantify 
the effect of selection and irreversible mutation respectively, the second term quantifies the effect of random 
drift. 

The stationary distribution can be expressed as 



Results 

Fixed points and their adaptiveness 

To understand the mechanism of Muller's ratchet, a full characterization of dynamical process is a pre- 
requisite for obtaining more accurate decaying time. Here we study the dynamical behaviors by investigating 
the position and adaptiveness of all fixed points. We further derive the parameter regions for all possible 
cases. 

According to general analysis of a dynamical system, letting 



$(x) = 



2Nfi(l - a) 

1 — CT/U 

| 2JV(1- M ) 

1 — (Jfl 



ln(l — x) — lnx(l — x) 



ln(l — a + xcr(l — ji)). 



(30) 




(31) 



$'(x) = 0, 



(32) 



we get 



2ct(1 - fi)(N - l)x 2 + (2N{p -a) + 3a-afj,- 2)x + (1 - a) = 0. 



(33) 



We solved the Eq.(32) and found two fixed points. If we denote 



a = 2 - 3a + an + 2Na - 2N >. 



(34) 



13 = 8a(l- M )(jV-l)(l-a). 



(35) 
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They are 



a T y/a 2 - (3 

X ^ = mhm ' (36) 



For two singular points x — 0, 1, if x — > 1, and cr G (/i, (2N(jl — l)/(2Nfj, — /i)), $(x) — > — oo. So the 
population is unadaptive at x = 1. When x — > 1, and cr e {{2Nfi — l/(2Nfi — fi), 1), $(x) -4- +oo. So the 
population is adaptive at x = 1. For a; — > 0, <&(x) — > +oo in almost parameters regimes except a = 1. So 
the population is always adaptive at a; = 0. When cr = 1, the viability of the suboptimal class is zero, so 
populations stay at the initial state, the corresponding minimum of adaptive landscape demonstrates the 
state with allele frequency x = 0. 

Here we address dynamical behavior by the positions of two real inequivalent fixed points x\ < X2 first. 

I) We find two different real fixed points in the regimes of fj, e (0,2/(2iV - 1 + 2^/N(N - 1))) and 
a E {/j,, 1); in the regimes of ^ e (2/(2JV - 1 + 2^/N(N - 1)), 1) and a <= ((2 + 2/j,-10Nn + 4N 2 fi + 2Nfi 2 + 
4(1 - fi) \J N(N - 1){{2N - - l))/(p - 2N + l) 2 , 1) except the regime of fi e {{2N - 1)/W(N - 1), 1) 
and cr = (2iV>- l)/(2iV>-/i). 

We discuss the position between them and the boundary points 0, 1 and adaptiveness of them in the 
following. 

i) 1 < xi < x 2 

In the regimes of fi e (1/(2N - 1),2/(2JV - 1 + 2^/N{N - 1))] and cr e (/z, (2JV/z - l)/(2N/i - /x)); in the 
regimes of ^ e (2/(2iV - 1 + 2^/N(N - 1)), (27V - l)/4iV(iV - 1)) and cr e ((2 + 2fj, - 10iV> + 47V 2 /i + 
2A^i 2 + 4(1 - n)y/N{N - l)((2N - 1> - l))/(^ - 2iV + l) 2 , (2W> - l)/(2iV/u - m)), the fixed points satisfy 
1 < xi < x 2 . At the same time the singular point x = 1 is unadaptive. There is one adaptive state with 
allele frequency x = in the system. Populations tend to evolve to the adaptive state. 

ii) 1 = xi < x 2 

In the regions of ^ e (1/(2N — 1), (2N — 1)/4N(N — 1)) and cr = {2Npt - l)/{2Npt - fi), the two fixed points 
satisfy x\ = 1, 1 < x 2 . The state with allele frequency x = 1 is unadaptive. There is one adaptive state with 
allele frequency x = in the system. 

iii) < xi < 1 < x 2 

In the regimes of \i e (0, 1/(2N — 1)) and a e (fi, 1); in the regimes of \i e (1/(2N — 1),1) and cr g 
((2Nfi — l)/(2JV/x — 1), the fixed points satisfy < xi < 1, 1 < x 2 . There is only one unadaptive state 
with allele frequency x = x\ in the system, and two adaptive states with allele frequency x = 1 and x = 
occur in the system. Populations tend to evolve to the adaptive state dependent on the position of the initial 
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state. If the initial state with allele frequency is greater than x\, populations tend to evolve to the adaptive 
state with allele frequency x = 1. 

iv) = xi < 1 < x 2 

In the regime of [i e (0, 1) and a = 1, the fixed points satisfy x\ = 0, 1 < a; 2 - When selection rate u = 1, 
the process lies at the initial state because for this case, the viability of the sub-fittest class is zero. 

v) — x\ < x 2 < 1 

The case = x\ < x 2 < 1 is impossible. For x\ = 0, that is a = 1, at the same time x 2 must be greater 
than one. 

vi) < x\ < x 2 < 1 

In the regimes of fj, e ((2AT-l)/4iV(iV-l), (2AT-l)/(4iV-3)) and a e {{2 + 2^-lON ^i+AN 2 ^i+2N ^i 2 
fJ.)y/N(N - 1)((2JV - 1)// - 1))/(m-2^+1) 2 , (2JV/i-l)/(2JV/i-^i)); in the regimes M e ((2/V-l)/(4/V-3), 1) 
and ct e ((2 + 2^ - 10JV/i + 4N*n + 2Nfi 2 + 4(1 - fi)y/N(N - 1)((2JV - \)fi - 1))/(m - 27V + l) 2 , (2iV,u - 
l)/(2Nfi — n)), the fixed points satisfy < x\ < x 2 < 1. The state with allele frequency x\ is unadaptive 
while that with allele frequency x 2 is adaptive. There are two adaptive states with allele frequency x = 
and x = x 2 and two unadaptive states with allele frequency x — 1 and x — x\ in the system. Populations 
evolve to which adaptive states dependent on the initial position. 

vii) x\ < or x 2 < 

The case x\ < is impossible, and the case x 2 < is impossible. 

II) Then we discuss the case of two equivalent real fixed points x 2 = x\. 

In the regimes of ^ e (2/(2iV - 1 + 2y/N(N - 1)), 1) and a = (2 + 2fi - lONfj, + 4N 2 ^i + 2N(i 2 + 4(1 - 
fj,)y/N(N- l)((2N - !))/(// - 2iV + l) 2 , we find two same fixed points 

Xl < 2 = mi-mV-i)- (37) 

i) 1 < X\ t 2 

In the regimes of fi e (2/(2N - 1 + 2^(^-1)), (27V - 1)/4N(N - 1)) and a = (2 + 2/z - 10/V/x + 47V 2 /i + 
2JV^ 2 +4(l- y u) v /JV(JV - 1)((27V- 1)^ - l))/(^-2iV + l) 2 , there are two same fixed points satisfy 1 < x lt2 , 
and they are unadaptive. There is one adaptive state with allele frequency x — in the process. 

ii) 1 = £1,2 

At the two points of ((27V - 1)/4N(N - 1), 2N/{2N - I) 2 ) and ((27V - l)/(4/V - 3), (4(iV - 1)(3 - 6JV + 
4iV 2 ) + 8(Ar-l)(4iV-3)^iV(iV- l)/(4Ar-3))/(2A-l) 2 ), there arc two same fixed points satisfy x ia = 1, 
and they are unadaptive. There is one adaptive state with allele frequency x = in the process. 
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iii) < x lf2 < 1 

In the regime of ^ e ((27V - 1)/47V(7V - 1), (27V - I) /{AN - 3)), a = (2 + 2^ - lOTV/x + 47V 2 + 27V/1 2 + 
4(1 - fJ,)y/NjN - 1){{2N - l)/i - 1))/(m - 27V + l) 2 ; in the regimes of ^ e ((27V - 1)/(47V - 3), 1) and 
a = (2 + 2/x - 107V^ + 47V 2 ^ + 27V^ 2 + 4(1 - \i)yjN{N - 1)((27V - l)/i - l))/(/x - 27V + l) 2 , there are two 
same fixed points satisfy < x\p. < 1, and they are unadaptive. There is one adaptive state with allele 
frequency x = in the process. 

Ill) Finally we consider two imaginary fixed points \x\ \ — \x 2 \- Where the |.| denotes the length for an 
imaginary points. 

In the regime of /x e (2/(27V - 1 + 2- v /7V(7V - 1)), 1) and a € (it, (2 + 2^ - lOTV^i + 47V 2 ^ + 2Nfi 2 + 4(1 - 
fi)y/N{N — 1)((27V — l)/x — 1))/(m — 27V + l) 2 ), there are two imaginary fixed points in the system. There 
is only one adaptive state with allele frequency x = 0. Populations always evolve to the adaptive state. 

Especially $'(a;) = is a linear equation, the fixed point is read out from the expression of f{x). We can 
measure the adaptiveness by the value of adaptive landscape. The bigger the value of adaptive landscape 
is, the more adaptive the corresponding state would be. The corresponding area of the fixed points in 
parameters plane is the Fig.l. 

Irreversible mutation, selection and random drift balance 

Concretely we divide mutation rates into three regimes. One is with mutation rates \i £ (0, 1/(27V — 1)) 
and selection rates a e (/U,l). Another is mutation rates [i e (1/(27V— 1), (27V— 1)/(47V(7V— 1))] and selection 
rates a € (1/(27V- 1), (27V> - 1)/(2j\> - //)). Another is with mutation rates fj, e ((27V - 1)/47V(7V - 1), 1) 
and selection rates a <G (27V/(27V — 1) 2 ,(2A^ — 1)/(27V^ — ix)). The first parameter regimes, a part of 
regions of I iii), corresponds to the case that mutation rates lie in the lower regime. And the adaptive 
landscape is U-shape. The second parameters regime, same of I i), is the middle regime. Adaptive landscape 
with these parameters demonstrates the monotonic decreasing. The third parameters region, same of I vi) , 
demonstrates the clicking process. The adaptive landscape of full parameter regimes is visualized as Fig. 2. 
From the expression and visualization of adaptive landscape $(a;), we may find there are two singular points 
and 1 of adaptive landscape, characterized by infinity values in Fig. 2. Singularity from the maximum of 
adaptive landscape indicates the population being adaptive while singularity from the minimum of adaptive 
landscape indicates the population being unadaptive. 

Fig. 2 demonstrates the whole process of the population evolution including the forming and losing the 
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fittest class except a = 1 demonstrated by red line in the middle graphs. Because the viability of the class 
with allele a is not zero and deleterious mutation is arbitrary. The adaptive state under the condition of 
the parameter a = 1 only means the initial states. With increasing selection rates the fittest class A forms 
quickly while with increasing mutation rates the fittest class A loses. In the lower mutation rates regime 
green line describes the population is likely to move to the fittest class with increasing selection rates, the 
process is dominated by selection. Blue line and red line manifest the losing process of allele A. Because the 
mutation rates are lower, selection rates dependent of mutation rates are lower, these factors result in the 
change of fittest class is not easy. In the end there are two adaptive states demonstrated by green line in the 
process. In the high mutation rates regime, blue line describes the population is likely to move to the fittest 
class so that the population exists in the form of coexistence of A and a. Because deleterious mutation rates 
are higher, as a consequence allele a occurs. Selection rates dependent of mutation rates tends to survive 
allele A. But the evolutionary process is dominated by the irreversible mutations, the fittest class A loses. 
There are two adaptive states in the process under the balance of irreversible mutation and selection. In the 
middle graph section, mutation rates and selection rates lie in the transition regime. It demonstrates the 
losing process demonstrated by red line. So we can draw the conclusion that the click process occurs when 
there are two stable states in the process. 



Characterization of the single click time 

We visualize the adaptive landscape, then one may wonder about how the population moves from one 
peak to another and how long it might be to move from one maximum to another. The process was first 
visualized by Wright in 1932. In addition, the problem of transition from metastable states is ubiquitous in 
almost all scientific areas. Most of previous works encounter finite potential barriers from the physical point 
of view. Interesting issue here is that we touch upon infinite potential barriers under the circumstance of 
well defined two adaptive states. Then we manifest the derivation of a single click time. The time of a click 



of the ratchet is recognized as the random time of loss of the fittest class 10 



The single click time is well defined when there are two fittest classes in the process. It means the 
interval time between extinction of the two fittest classes. The corresponding processes are that there are 
two well-defined adaptive states in the system. Corresponding graphs of adaptive landscape is Fig. 3. 

To evaluate the single click time and show the further power of adaptive landscape, in the following we will 
demonstrate how the single click time from one adaptive state to another is derived in this framework. After 
straightforward calculation, backward Fokker- Planck equation corresponding to Eq.(27) can be expressed 
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with the property of time homogeneous in the following form 29 30 

d tP (x, t) = (/(a;) + eD'(x))d x p(x, t) + eD(x)d 2 xP (x, t). (38) 
General single click time dependent on initial Dirac function satisfies 

(f(x) + eD'{x))d x T{x) + eD{x)d 2 x T{x) = -1. (39) 

With 

T'(1) = 0, 

T(0) = 0. (40) 

Above treatment is valid. Because populations evolution is according to Muller's ratchet, that is in the 
presence of deleterious mutation, without any recombination, but with selection and random drift. And the 
model in discrete manner demonstrates the transition probabilities are from the boundary x — to its 
next. So the boundary x = only absorbs flux from its next, the boundary is absorbing. The probabilities 
from boundary x = 1 is not zero because /j ^ and /i, er ^ 1. And the population can not be out of the 
boundary x = 1. So the boundary x = 1 is reflecting. The general solution corresponding to Eq.(39) is 

T{x) = J ^-rexp(-$(y))dy j* ' exp(*(«))dz, (41) 

here ${x) = f \f{x')/ 'D{x'))dx' '(e - 1). 

Here the evolutionary process occurs when x € [0,1]. We are more interested in the transition time 
between the two adaptive states x = and x = 1. In the process, there are two important states a;*, £g. 
Interval (0, 1) contains a potential well at x* and a potential barrier at Xq. The single click time is composed 
of two elements, one denotes forming process of fittest class, the other describes losing process of fittest class. 
In general, the time spent on forming process is much smaller than that spent on losing process. So the 
transition time approximates to the time spent on losing process. Because we assume that near Xq we can 
write 

*(a:) ^(x*)-^^) 2 . (42) 

and near x* 

4>(z)« $(z*) + i(^>. (43) 
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At the same time, if the central maximum of $(x) is large compared with 1/N, then exp($(z)) is sharply 
peaked at Xq, while exp(— $(y))/D(y) is very small near y = x* . Eq.(41) is evaluated as 

Ti^o ~ J j^-j exp(-$(y))dy J exp($(z))dz 
27ra / ^'exp($(^) - $(x*)) 

^i^exp^S) -$(**)). (44) 



oc 



From the expression of Eq.(44), the single click time is not sensitive to the assumption of Eq.(40). 

In the higher mutation rates regime, where Xq approximates to an adaptive state which is near enough 
to 1, x* corresponds to the unadaptive state that the population lies between the adaptive states and Xq. 
The potential barrier is finite. According to classical derivation corresponding to Eq.(44) the single click 
time approximates to 



' 1— s-0 



= lim / -F^ cx P(-®(y))dy [ cxp($(z))dz 
x^iJo eD{y) J y 

f x " (1 _ y )2iV M ((T-l)/(l-<r/i) 

~ 2N Jo (1 - o + ya(l - M ))2W(i-ri/(i-M dy 



fz-\l- 



\((TAi-l-2A^((T-l))/(l-(Tp) 



x (1 - a + az(l - fiffW-MV-'ridz 

^L- )CM ^ x *)-<P(x*j) 

N(N — l) 2 cr 2 (l — fi) 2 



(45) 



(a- V" 2 - P)(4*( N - 1)(1 - M) - " + V a2 - P) ' 
Here Xq is the fixed point X2, and x* is the fixed point x\, parameters a, (3 is the same as Eqs.(34) (35) 
respectively. The difference of potential is 

*(*S) -*(**) = 2iV/ ^±^ lnfl-^Vlnfl + ^ 

\ — a\i \ 1 — xi ) \ x\ 

, 2iV(l - M) ln L | <7(l-M)(x 2 -^i) 
1 — <j/U \ 1 — (J + x\o(\ — /i) 

2JV M (1 - <r) ^ / _ 2^^^ ^ 



ln | 1 



4(a 2 - /?) 



(a - y/a 2 - /3) 2 J 

+ 2N(1^)J 1+ \ 

I -vfi \ AN -6- 6Na + 7cr — <T/U + 2iV/i + V" 2 ~ /? / 
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Where a and f3 are the same as Eqs.(34) and (35)The approximated single click time varies with mutation 
rates in Fig. 4. The single click time Ti^o increases with population size N in certain regime, decreases 
with mutation rates fi and selection rates a in the parameters regime [i e (2N/AN(N — 1),1) and a e 



((2 + 2fi-10N^i + 4:N 2 ^ + 2Nfi 2 +4:(l-fi)^/N{N - l)((2N - l)fi- l))/{p-2N + l) 2 , (2JV>-l)/(2iV>-/i)). 



Because in the regime, with selection rates increasing, the difference of potential between two fixed points 
decreases as the demonstration in Fig. 2, the viability of suboptimal class decreases, populations evolve to the 
fittest class. These results in the the single click time shorter. In another hand, with deleterious mutation 
increasing, the population of suboptimal class increases, the difference of potential between two fixed points 
decreases as the demonstration in Fig. 2. These also results in the single click time shorter. 

For the lower mutation rates regime, where the potential barrier is infinite. The single click time can be 
estimated also, x* corresponds to the fixed point x\ that the population lies at the lowest potential. 



From expression of Eq.(47), the single click time goes to infinity with mutation rates tends to zero in the 
parameters regimes of \i e (0, l/(2iV— 1)) and a g (fj,, 1). And when parameters regions lie \i € (1/(2N— 1), 1) 
and a g ((2N/J, — 1)/(2N(/, — fi), 1), the results of the single click time is not sensitive to the population 
size. Biologically if deleterious mutation accumulates, the viability of suboptimal class increases as the 
demonstration in Fig. 2, these results in the single click time longer. 




x (1 - a + az{l - ii)) 2N ^-^l {1 -°^dz 



1 — 071 



(47) 
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Analogous to the derivation of Ti^ , we can calculate 



( x 1 r 

T ^i= lim / — — exp(-$(y))dy / exp($(z))dz 

/■x _ ^2iV M (<r-l)/(l-<vO 

~1 1 ™0 2N J 1 (1 - a + - a1 ))2AT(i- m )/(i-<t /1 ) d f 



(cr^-l-2JV J u(cr-l))/(l-CT ( u) 



x / ^(l-z 
x (1 - a + <rz(l - /1 )) 2JV ( 1 -")/(i-"M) ( fe 

~ 7 (1 - a + - M ))2iV(i- M )/(i-^) a V 
rV 

x / + . . .)dz 

Jo 

— > oo. (48) 

Compared with the singular point a; = 1, the difference between two singular points x = and .t = 1 is 
the mutation rates /u. The mutation rates /U > at x = 1, this results in the power of (1 — z) is not negative, 
so the single click time is finite. The infinity of the single click time from x = comes from that the mutation 
from unfavored allele a to favored allele A is zero, this results in the second integral nonintegrable because 
of the negative power of z. These results are consistent with biological inference. Biologically because the 
absence of back mutation, the accumulation of deleterious mutation, once the population arrives at the state 
that almost individuals are with allele a, the population is absorbed to the state and can not leave with 
high probability. Mathematically, because the second integral is singular for the singular point x = 0. And 
the integrated function is a fraction respect to argument x, but the highest power of denominator is smaller 
than that of numerator. That results in the power of x is 2N — 2. As a consequence the second integral is 
singular. 

Discussion 

We analytically construct adaptive landscape. The constructive method is independent on the existence 
and normalization of stationary distribution. We demonstrate the position and adaptiveness of all fixed 
points for the whole parameters regimes under the condition of the diffusion approximation. An interesting 
thing is the imaginary fixed points occurring. We give the parameters regions of their occurrence. However, 
we have not found any study of Muller ratchet for the fixed points to give a complete description. In addition, 
we give the description of escape from infinite potential. However, intuitively infinite potential means the 
population lies at adaptive state. The transition from the adaptive state can not occur. Here we find that 
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the escape from infinite potential can not occur when the boundary is absorbing. So we define the absorbing 
boundary by adaptive landscape and the single click time without any extraneous assumptions. 

The model with discrete manner describes the nature of populations evolution. Here we give two special 
cases. One is that the population lies at that state with allele frequency x = 0, the other is that the 
population lies at the state with allele frequency x = 1. We compare the model with discrete and continuous 
manners to conclude the definition of boundary conditions without any extraneous assumptions. The model 
with continuous manner is derived under the condition of enough small space change Sx in one generation, in 
another word, when population size N is much bigger. The model with continuous manner can correspond 
to the model with discrete manner. 

When allele frequency x = 0, the matrix for transition probabilities demonstrates that the population 
only absorbs the flux from the next, and it can not output any flux. In the model with continuous manner 
the value of adaptive landscape with allele frequency x = arrives at the maximum except a = 1. This 
demonstrates the state with allele frequency x = is almost always adaptive. Then we calculate the single 
click time from the boundary with allele frequency x — 0. We find the single click time is infinite. So we 
draw the conclusion that the state with allele frequency x — absorbs flux, and does not output flux. The 
boundary is absorbing. When the allele frequency x = 1 in the model with discrete manner, the transition 
probability is 

(iV !( M (1 - a)) N - m (l M) m )/H(iV - m)!(l - *(i) N ), 

that is, when /i > 0, the transition probabilities are not zero, the state with allele frequency x = 1 can 
input any flux to its next state. The single click time from the state is finite in the model with diffusion 
manner, however the potential at the state could be infinite, then we draw the conclusion this boundary is 
not absorbing. This is consistent with the biological understanding. 

This article presents an approach to estimate the single click time of Muller's ratchet. Furthermore, it 
define the absorbing phenomenon by the single click time without any extraneous assumptions. Inspired 



by 16 10 , we connect Muller's ratchet to one locus Wright-Fisher model with asexual population including 
N haploid individuals. And our model is represented as a Fokker-Planck equation. We give a complete 
description for the position and adaptiveness of all fixed points in the whole parameters regimes. This 
is first done bases on diffusion approximation. The investigated elements are at the allele level. This 



is different from Ref. 10 . Our method does not need the existence and normalization of the stationary 
distribution. Our constructive method is independent of the stationary distribution. Compared with the 
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method based on diffusion approximation [15] |2j, mathematically it is described as stochastic differential 
equations. Our method investigates the global dynamical property of the system, and reduces the complexity 
of calculating stochastic differential equations. In addition, the boundary condition of these stochastic 
differential equations is prescribed. Compared with Ref. [To], They added Dirac function to the boundary. 
But this is not appropriate for the adding non-differential Dirac function to stationary distribution, and 
stationary distribution should satisfy diffusion equation. However, the treatment is convenient for computing 
the stationary distribution. The stationary distribution of theirs is equivalent to our adaptive landscape. 
They had not given the shape of adaptive landscape when the mutation rate lies in the lower regime. We 
use the model defined in the interval (0,1) to describe the absorbing boundary. We check the biological 
phenomenon by the model with both discrete and continuous manners. This is a new method to handle the 
boundary condition. We investigate the absorbing phenomenon by it without any extraneous assumptions. 

To summarize, we have obtained two main sets of results in the present work. Most importantly, we 
find that the absorbing phenomenon is characterized by the adaptive landscape and the single click time 
without any extraneous assumptions. First, we demonstrate the adaptive landscape can be explicitly read 
out as a potential function from general diffusion equation. This not only allows computing the single click 
time of Muller's ratchet straightforward, but also characterizes the whole picture of the ratchet mechanism. 
The adaptive landscape has rich structures such as finite and infinite potential, real and imaginary fixed 
points. We analytically demonstrate the adaptive and unadaptive states for the whole parameters regimes. 
We find corresponding parameters regimes for different shapes of adaptive landscape. Second, we give the 
formula about the single click time with finite and infinite potential. And we find the single click time with 
finite potential increases with selection rates and population size increasing, decreases with mutation rates 
increasing, the single click time with infinite potential is insensitive to the population size. These results give 
a new understanding of infinite potential and allow us a new way to handle the absorbing phenomenon. In 
this perspective our work may be a starting point for estimating the click time for Muller's ratchet in more 
general situations and for describing the boundary condition. Such demonstration suggests that adaptive 
landscape may be applicable to other levels of systems biology. 
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mutation rate n 

Figure 1: Relation of fixed points and parameters for the system in all regimes. The regime represented by III i) with 
parameters regions a 6 (<7i, 02) and /x G (2/(2N — 1+2-y/ N(N — 1)), 1) has two imaginary fixed points. The red curve 
denoted by II) with parameters satisfying a = 02 and fi G (2/{2N — 1 + 2sJ N(N — 1)), 1) has two equivalent fixed 
points. The three cases of II) occur in the intervals. The regimes denoted by I) with parameters satisfying a G (02, 0-5) 
and n G (2/(2N - 1 + 2^(^-1)), (2N - 1)/AN(N - 1)); satisfying o G (a 2 , o~ 5 ) and fj, G ((2N - 1)/AN(N - 1), 1); 
satisfying a G (0-1,0-5) and fi G (0, l/(2iV-l)); satisfying a G (0-1,0-5) and ^ G (1/(2N- 1), 2/(2N- 1 + 2 y/N(N - 1))) 
have two real fixed points. The five cases of I) occur in the regions. 
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mutation rate \i 

Figure 2: Adaptive landscape against allele frequency x with mutation rates and selection rates in all regimes, x 
label represents allele frequency, fi label represents mutation rates while vertical label is corresponding value of $. 
Assume population size is constant 7V=100. In the right graph section, green line corresponds to the case of I iii), 
where parameters satisfy fi=0. 00005, <r=0.0005. Blue line approximates to the case of III, where parameters satisfy 
/i=l/198, 0"=1/198. Red line corresponds to the case of I i), where parameters satisfy /i=0. 0050252, a~0. 00505. In the 
middle graph section, green line represents the case of I ii), where parameters satisfy fi=0. 0050252, a~0. 005039925. 
Blue line represents the case of II ii), where parameters satisfy /i=0. 0050252, <r=0. 0050503775. Red line represents 
the case of II i), where parameters satisfy /i=0. 005025, <r=0. 00503. In the left graph section, green line represents 
the case of, where parameters satisfy /i=0.04, <r=0.06 II iii). Blue line represents the case of I vi), where parameters 
satisfy )i=0.04, <r=0.19. Red line represents the case of I iv), where parameters satisfy /i=0.04, cr=0.99. 
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mutation rate |i 



allele frequency x 



Figure 3: Adaptive landscape with two adaptive states. iV=100, blue line represents /i=0.04,cr=0.19 corresponding 
the case for I vi) while green line stands for /j, = 0.000005, a~0. 00005 corresponding the case for I iii). 
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■ high mutation regime, u=21 0/(400*99) 

■ low mutation regime, u=21 0/(400*99) 

■ high mutation regime, u.=500/(400*99) 
low mutation regime, u=500/(400*99) 
high mutation regime, u=1 000/(400*99) 
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The approximated single click time decreases with mutation rates increasing and selection rates 



Figure 4: 
decreasing in the regime 
fj,)y/N(N-l)((2N 



jj_ e {2N/AN{N - 1),1) 
1))/(M- 2N + l) 2 ,(2iV M 



and a 6 
- l)/(2Nn - 



((2 + 2^ - WNfi + 4A?> + 2Nfj? + 4(1 



1))/(M - 27V + iy,(2N(j, - 1)/(2N(jl - fi)) denoted by I iv). Assume 7V = 100, blue, 
black and dashed lines, corresponding high mutation rates regime, represent fi = 210/(400 * 99), n = 500/(400 * 99) 
and /x = 1000/(400 * 99) respectively. Red and green lines, corresponding low mutation rates regime, represent 
H — 210/(400 * 99) and n = 500/(400 * 99) respectively. The approximated single click time increases with mutation 
rates increasing in the regime /i g (l/(2iV — 1), 1) and a £ ((2Nfi — l)/(2Nfi — /x), 1) denoted by I iii). Green dotlinc 
represents this case. 
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